Interaction-induced charge and spin pumping through a quantum dot at finite bias 
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We investigate charge and spin transport through an adiabatically driven, strongly interacting 
quantum dot weakly coupled to two metallic contacts with finite bias voltage. Within a kinetic 
equation approach, we identify coefficients of response to the time-dependent external driving and 
relate these to the concepts of charge and spin emissivities previously discussed within the time- 
dependent scattering matrix approach. Expressed in terms of auxiliary vector fields, the response 
coefficients allow for a straightforward analysis of recently predicted interaction-induced pumping 
under periodic modulation of the gate and bias voltage [Phys. Rev. Lett. 104, 226803 (2010)]. We 
perform a detailed study of this effect and the related adiabatic Coulomb blockade spectroscopy, 
and, in particular, extend it to spin pumping. Analytic formulas for the pumped charge and spin in 
the regimes of small and large driving amplitude are provided for arbitrary bias. In the absence of 
a magnetic field, we obtain a striking, simple relation between the pumped charge at zero bias and 
at bias equal to the Coulomb charging energy. At finite magnetic field, there is a possibility to have 
interaction-induced pure spin pumping at this finite bias value, and generally, additional features 
appear in the pumped charge. For large-amplitude adiabatic driving, the magnitude of both the 
pumped charge and spin at the various resonances saturate at values which are independent of the 
specific shape of the pumping cycle. Each of these values provide an independent, quantitative 
measurement of the junction asymmetry. 

PACS numbers: 72.25.-b, 73.23.Hk, 73.63.Kv 



I. INTRODUCTION 

The generation of a dc current through a mesoscopic 
system is usually associated to a bias voltage maintained 
between the contacts to the exterior world. Remark- 
ably, charge and spin transport can even be achieved 
in the absence of an external bias by the cyclic modu- 
lation of some of the parameters of the system^ When 
this modulation is slow compared to the characteristic 
dwell time of the electrons, the transport mechanism 
is called adiabatic pumping. Here, the pumped charge 
is of geometric nature, since it depends on the specific 
shape of the path sustained by the system's parame- 
ters but not on its detailed time evolution^— For ap- 
propriate modulation setups^ii the pumped charge af- 
ter one period may be quantized in units of the elec- 
tron charge^— motivating its use as a highly precise 
current standard for quantum metrology-^ or in the 
initialization^ and processing of coherent states^ 3 - - — in 
the context of quantum information. In the opposite 
limit where the transferred charge is not necessarily quan- 
tized, the pumping mechanism is dominated by quan- 
tum interference^ of the coherent electrons in the device. 
In the last years, adiabatic pumping was widely studied 
both experimentally-^— and theoretically. In noninter- 
acting systems, a well-established theory was formulated 
by Brouweri^i It makes use of the concept of emissiv- 
ity which was introduced in the scattering matrix ap- 
proach for time-dependent systems at low frequency by 
Buttiker, Thomas and Pretre. 22 Importantly, this formal- 



ism is adequate as long as interactions can be described 
on a self-consistent mean- field level^ Within this formal- 
ism, several aspects of adiabatic pumping were explored, 
covering diverse effects such as dissipation and noise 
or spin polarized pumping^ Further works dealt with 
different setups including normal metal-superconducting 
heterostructures^ pumping by surface acoustic waves^ 
and graphene-based quantum pumps.— 

Pumping through confined electron systems domi- 
nated by a strong Coulomb interaction is a particularly 
challenging topic since the mean-field approach breaks 
down and a new formulation is necessary. Several stud- 
ies addressed interaction effects in specific setups and 
regimes In Ref. [ID, pumping is investigated in inter- 
acting quantum wires. By using a slave boson mean field 
approximation, Aono studied adiabatic pumping through 
a quantum dot in the Kondo regime* 3 ^ This regime was 
also treated in the Toulouse limit^i and for nonadia- 
batic pumps^ 2 . by using the Keldysh Green's function 
technique. Pumping through open quantum dots was 
described by employing bosonization techniquesj 33 i In 
the Coulomb blockade regime, spin pumping was ad- 
dressed through a numerical calculation of the reduced 
density matrix of a double quantum doti 3 ^ An expres- 
sion for the adiabatic pumping current in interacting sys- 
tems was derived usin g a noneq uilibrium Green's func- 
tion technique in Refs. |36| andl37l A diagrammatic real- 
time approach^ was used to investigate several aspects 
of adiabatic pumping through weakly-coupled interact- 
ing quantum-dot systems^— and served as the basis for 
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non-equilibrium renormalization group studies that treat 
the tunneling non-perturbatively4£ 

Among the above mentioned studies, only a few dis- 
cuss the modulation of the applied bias4£ In particular, 
pumping around a nonequilibrium working point induced 
by a static nonlinear bias was addressed in noninteracting 
systems.— Recently, a strongly interacting single-level 
quantum dot with a modulation of the gate and bias 
voltage was investigated^ For this modulation setup, on 
top of a dc current produced by the bias, an additional 
adiabatic dc current is generated by the Coulomb inter- 
action. Interestingly, this interaction-induced pumping 
current can be accessed by using lock-in techniques, and 
was suggested as a new spectroscopic tool to probe in- 
ternal properties of the system, like spin degeneracy and 
junction asymmetries. Similar effects were reported^ for 
an open quantum system when controlling the tempera- 
tures and chemical potentials of the reservoirs. In Ref.0, 
the zero-frequency pumping noise in adiabatically driven 
quantum dots is discussed for time-dependent bias, re- 
vealing further information on the tunnel coupling asym- 
metry in cases where the pumped charge is zero. 

In this paper, we investigate interaction-induced 
charge pumping in detail and extend it to the spin de- 
gree of freedom. We focus on the interplay between 
the strong local Coulomb interaction in a quantum dot 
and the nonequilibrium effects induced by finite bias and 
the modulation around this working point. To describe 
the dynamics of the local system, the coupling to the 
leads and the frequency of the modulation are treated 
perturbatively.^ 8 - In particular, we restrict ourselves to 
the single-electron tunneling (SET) regime. 

We derive a general expression for the adiabatic charge 
and spin currents in response to a change in the internal 
occupations of the dot induced by the driving param- 
eters. The coefficients of this response of the current 
are related to the emissivities to the leadsj^ establish- 
ing a connection to Brouwer's pumping theory 21 which 
is applicable as long as the single-particle picture holds. 
Moreover, we write the pumped charge and spin as the 
flux generated by auxiliary vector fields in the space of 
the parameters. This allows us to find the general con- 
ditions under which a finite pumped charge or spin may 
occur, independent of details specific to the model or the 
pumping cycle. 

We apply this strategy for the case of a modulation 
of the gate and bias voltage and analyze the influence of 
the local interaction on the generation of pumping. The 
recently introduced "stability diagram" for the pumped 
charge^ is shown to be readily understood in terms of 
these vector fields and extended to pumped spin. A de- 
tailed analysis of these diagrams is provided for all the 
regimes of the applied bias, including analytic fitting for- 
mulas for the pumped charge and spin in the limits of 
weak and large driving amplitude. We find that, al- 
though not quantized, the pumped charge (and spin) 
for large driving amplitudes saturates at plateau values 
which depend on the asymmetry in the coupling to the 
leads and the working region of the applied voltages. At 



zero magnetic field, a surprising and simple relation be- 
tween the pumped charge at high-bias and at zero-bias 
is found. It is exploited for a quantitative determination 
of the junction asymmetry by two single measurements. 
In a finite magnetic field, the spin and charge pumping 
are not anymore trivially connected at high-bias. In this 
regime, pure spin pumping occurs when the coupling to 
the leads is symmetric. 

Our paper is organized as follows. In the next sec- 
tion we introduce the model and review the theoretical 
framework used in the calculation of the pumped charge 
and spin. In Sec. IIII1 we consider the pumped charge 
and spin for the specific modulation of the gate and bias 
voltage, first concentrating on the role of the interaction 
and then, in Sec. IIIID| we discuss the effects induced by 
an external magnetic field. We summarize our results in 
Sec. ED 

II. MODEL AND FORMALISM 
A. Model 

We consider a single-level quantum dot with Coulomb 
interaction weakly coupled to two noninteracting leads 
as sketched in Fig. [TJ The full system, containing the 
dot, the left (L) and right (i?) leads, and the tunneling 
between dot and leads, is described by the total Hamil- 
tonian H(t) = iTdot(^) + H les (t) + H tun . The quantum 
dot Hamiltonian is given by 

Hdotjt) = e a {t)h a + Uhfhi, (1) 

where we denote the spin-resolved number operator by 
n a = d\d„ and U is the Coulomb charging energy. The 
fermionic operator d\ (d a ) creates (annihilates) an elec- 
tron in the dot with spin a =f, \.. The many-body eigen- 
states of the dot are characterized by their charge num- 
ber and spin by |0) for an empty dot, \a) — d\\Q) for a 
singly occupied dot with spin a and |2) = d]dt |0) for a 
doubly occupied dot. Their energies are then 0, e a and 
e t + e J. + ^' respectively, where e a (t) — e(i) — aB/2. The 
level position e(t) = —aV g (t) is capacitively modulated 
by the time-dependent gate voltage V g (t), with lever arm 
a < 1. Furthermore, B accounts for the Zeeman split- 
ting produced by an external magnetic field B = Be z 
in units e = h = fcs = gfis = I- Here we use a = ±1 
as a convenient notation for spin f and I, respectively. 
The underlying capacitive description of the parameters 
U and a is sketched in Fig.[2]and will be discussed below. 

The leads are described as reservoirs of noninteracting 
electrons through the Hamiltonian 

H les (t) = ^[e r fc + Hr(t)}cl kcr c rka , (2) 

rka 

where c\ k(J (c r k a ) creates (annihilates) an electron in the 
lead r = L,R with spin a J, and state index k. The 
eigenenergies of the leads are uniformly shifted by the 
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FIG. 1. Scheme of the considered model. The interacting 
quantum dot (grey circle) is capacitively coupled to two leads 
via Tl and Tr. The transport is controlled through time- 
dependent gate and bias voltages. The lines in the dot in- 
dicate the transition energies ejv — ejv_i, N — 1,2, for zero 
external magnetic field B. 



time-dependent bias voltage Vb(t) such that the electro- 
chemical potentials read fi r (i) = ±T4(t)/2 for r = L,R. 
These reservoirs are furthermore characterized by a tem- 
perature T. 

Finally, the tunnel coupling between the dot and the 
leads is determined by the tunnel Hamiltonian 

fftun = / ] Vrdjcrktr + h.C, (3) 
rka 

with the tunnel matrix element V r , which we assume to 
be independent of k and a. The tunnel-coupling strength 
T r = 2ir\V r \ 2 v r characterizes the rate at which tunnel 
processes take place. Here v r is the density of states in 
the r-lead, which is assumed to be energy-independent 
and with a band cutoff D r , which is the largest energy 
scale. 

We are interested in the simultaneous modulation of 
gate and bias voltages. These are driven around a time- 
independent working point specified by V g and V& with a 
fixed relative phase Acj> = <j> g — fax 

V x (t) = V X + SV X sm(nt + (p x ), x = g,b, (4) 

where SV X is the driving amplitude. At time to, we 
"switch on" the coupling between the dot and the leads 
and calculate the time-dependent nonequilibrium steady 
state at a much later time t. This includes two sources 
of nonequilibrium: the finite bias and the adiabatic driv- 
ing. On top of the stationary current flow, an additional 
time-dependent current, which can have a dc component, 
is generated by the periodic modulation. We work in the 
adiabatic regime where the driving period T = 2tt /VL is 
larger than the typical time spent by the electrons in- 
side the quantum dot. Here, the modulation frequency 
and the driving amplitude are limited by the adiabaticity 
condition aSV g ,SVb <C T /fi. 



1. Observables 

In the description of the dot occupancies through the 
reduced density matrix, we assume the reservoirs to be 
always in equilibrium. Since leads and dot are decou- 
pled at the initialization time to, the total density matrix 
is factorized in both subsystems as p — p rcs p, with the 
density matrix of the leads given by the one of a grand 
canonical ensemble 

i5 re8 = nfe-[^-MO^]/T (5a) 

r r 

Zr = Tre-l H ^-^ t '^ T , (5b) 

r 

and t' < to. Here N r — J2ka c rka Crkcr ^ s ^ ne P ar ticle- 
number operator in the r-lead and the temperature T 
is assumed to be the same in the two leads. Notice 
that the density operator of the reservoirs remains time- 
independent, meaning that the occupation in the leads is 
not affected by the modulation of the bias voltage. For 
electric field modulations with a frequency well below 
the plasma frequency of the leads (tens of THz in typical 
doped semiconductors) this assumption is well justified. 50 
The time evolution of the expectation value of an ar- 
bitrary operator R is formally obtained by 

R(t) = (R)(t) = TrTr(Rp(t)). (6) 

dot res V / 

In the next sections, we will describe the tunneling cur- 
rents I r (t) = (I r ) (t) and J r (t) = (J r ) (t) related to the 
charge and the spin component along the external field 
B, respectively, entering the r-lead. Since in the uncou- 
pled system Ho = -ffdot + H ICS the number of particles 
and the z-component of spin are conserved, the operators 
related to these observables are given by 

i r =i[H tun ,N r ], (7a) 
J r = i[H tun ,S r z ], (7b) 

with §1 — J2ka l c lka Crka '- Since we consider a system 
that preserves rotation symmetry around the magnetic 
field (z-) axis, the x and y components of the spin vector 
may only occur as a transient effect and are not required 
here. 



B. Real-time diagrammatic approach 

In this section we outline the theoretical framework 
used to identify the adiabatic contribution to the time- 
resolved currents of Eq. ©. As we will show in Eqs. (fT4")) 
and (I15|) , the adiabatic current can be interpreted as the 
delayed response of the dot occupation probabilities to 
the driving of the external parameters. The full rele- 
vance of this result will become clear in the next section. 
Eq. (|15p allows to apply the concept of emissivity to a 
system with strong Coulomb interaction and arbitrary 
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bias. Originally^ the emissivity was introduced to de- 
scribe capacitive effects on time-dependent response on 
a self-consistent mean-field level and was further used on 
adiabatic pumping through non- or weakly interacting 
systemsi 21 ' 37 

We start with the description of the relevant part 
of the dot's reduced density matrix, namely, its diag- 
onal elements, obtained after tracing out the degrees 
of freedom of the leads. The time evolution of the 
dot occupation probabilities, represented by the vector 
P(*) = (Po(t),P^{t),Pi(t),P2(t)) T , is governed by the gen- 
eralized master equation^ 



dt 



P(t) 



dt'W(t,t')p(t'). 



(8) 



Here, the change in the dot occupation probabilities, due 
to electron tunnel processes between dot and leads, is 
accounted by the kernel W(i, t'). In terms of the real- 
time diagrammatic technique developed in Ref. l5ll this 
kernel collects all irreducible diagrams in the Keldysh 
double contour. Its matrix elements W miTl (t,t') describe 
the transition from a state |n) at time t' to a state \m) 
at time t. It is important to notice that the transport 
properties are completely determined by the diagonal el- 
ements of the reduced density operator. In the chosen 
basis the off-diagonal elements, related to coherent super- 
position of different states, are decoupled from diagonal 
ones due to charge and spin conservation in the tunneling 
and therefore do not affect the currents. 

Since we consider a weakly coupled system bounded to 
a slow modulation of the energy levels, it is sufficient to 
describe the lowest order contribution in the tunnel cou- 
pling and in the time-dependent perturbation introduced 
by the driving. The occupation probabilities are thus ex- 
panded in powers of fl by p(t) = pj + Pt Q \ bearing in 
mind that fl/T <C 1. The first term (zeroth-order in fi) 
represents the instantaneous occupations and describes 
the steady state solution when the parameters are frozen 
at time t. Here the index t indicates the parametric time- 
dependence through the driving parameters {x(t)}, i-e. 

p^ — P^({x(*)})- The instantaneous occupations are 
obtained from the time-dependent kinetic equation in the 
stationary limit 



= W 



(<)_(*) 
Pt 



(9) 



together with the normalization condition e^p^ = 1, 
where e = (1,1,1,1) T , and we introduced the zero- 
frequency Laplace transform of the instantaneous kernel 
W t W = dt'W^it - if). In the SET regime we con- 
sider here, characterized by the linear dependence on T 

(i) 

of W; , the result coincides with the one from Fermi's 
golden rule. The next-to-leading term (linear in Q), p| a \ 
obeys the adiabatic correction to the kinetic equation 

(10) 



d (i) 
-Pt 



dt 



W(i) (a) 
1 Pt ■ 



The retardation correction p^ is determined by the 
competition between the driving (left hand side) and 



the inverse response times contained in (right hand 
side). The dot occupation probabilities are obtained by 
solving Eqs. ^ and (|10l) together with the normalization 

condition e T p^ = 1 and e T p[ a ^ — 0. From Eq. (fTUj) . the 
adiabatic corrections to the occupation probabilities are 
written in terms of the instantaneous contributions by 



» 



W 



d (i) 
-Pt > 



dt 



where the (invertible) matrix 

















ij 




ij 




ii 



(11) 



(12) 



includes the normalization condition e T pj = 0. 

The charge and spin currents in Eq. © need to be 
equally expanded in both the frequency and the tunnel- 
coupling strength V. The resulting observables are then 
split into instantaneous and adiabatic correction terms 



R 



(i) 

where W R t is the instantaneous kernel of the corre- 
sponding current R which, in the present approximation, 
is linear in T. We describe R^ by a scalar product with 
the time-derivative of the dot state occupation probabil- 
ities 



(i/a) 



(R) 



(i/a) 



R,i 



Pt 



(13) 



(i) 
R,; 

t d 



W 



-i d 



dt Pt 



(14a) 
(14b) 
(14c) 



with the sum running over the dot eigenstates, i.e. j = 
0,t,i, 2. Applied to the adiabatic charge and spin cur- 
rents, this equation defines the adiabatic current R as 
the response to a time-dependent variation in the instan- 
taneous occupation probabilities induced by the external 
modulation. The response coefficient, 



or: 



dp 



(15) 



determines the ratio at which the current R flows into 
the r-lead due to a variation in the occupation of the 
state j. The relevance of these coefficients lies in the 
fact that they distinguish the amount of charge (or spin) 
that enters into each one of the leads. As compared to 
the instantaneous solution, the response coefficients give 
information about the characteristic delay time for the 
current R£& 



C. Adiabatically pumped charge and spin 

Now that we have an explicit expression for the adi- 
abatic current, we can determine the charge and spin 
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pumped through the quantum dot during one modula- 
tion cycle. The purpose of this section is to relate the 
resulting pumped charge and spin to the emissivity of the 
contacts^, a well-known concept from scattering theory. 
It measures the amount of charge entering the r-lead due 
to the variation S\ of the driving parameter x- Following 
the reasoning by Biittiker et al^ for a slow variation of 
X, the charge entering the r-lead is related to the emis- 
sivity dN(r)/dx by 



oQi r = — -, — ox- 
ax 



(16) 



We are interested in the adiabatic 52 pumped charge when 
varying two parameters Xi(^) an d X2W over a cycle of 
the driving. Calculated as the integral of Eq. (fT6]). this 
reads^i 



T 



fdN(r)dxi dN(r)dxz 



dt . 

V d Xi 



dt 



dx2 dt 



(17) 



In terms of the adiabatic current of Eq. (|14p. this last 
can also be written as 



T 



dtl 



(a) 
r,t 



(18a) 
(18b) 



Since Pj t depends on t through the driving parameters, 
we rewrite its time-derivative in terms of Xi.2- A compar- 
ison with Eq. (|17|) allows one to relate the emissivity with 
the response coefficients and the occupation probabilities 
by 



dN{r) 
dx 



E 



dx 



(19) 



where x is either xi ° r X2 and the index t is removed to 

J (i) 

emphasize that ip? and p, are functions of x rather than 
t. Analogously, the above relation can also be generalized 
for the spin emissivity^ dS(r)/dx in terms of the spin 
current response coefficients as follows 



dS(r) 



E 



dp 



(i) 



dx 



(20) 



Written in this way, the emissivity is the weighted rate 
of change in the occupation probabilities due to the ex- 
ternal perturbation. The response coefficients tp 1 ^ and 

ifj r describe the rate at which the charge and the spin, 
respectively, are transferred to the leads when the occu- 
pation probabilities are changed by the driving parame- 
ters. The above Eqs. (fT9|) and (|20|) . also shown by Sela 
and Oreg2£ for adiabatic transport at equilibrium, extend 
the known result from scattering matrix theor y 21 ' 22 to a 
system with strong Coulomb interaction driven around a 
nonequilibrium steady state. 



1. Pseudo vector potential and pseudo magnetic field 

We now describe the pumped charge and spin in terms 
of auxiliary vector fields defined in the space of the driv- 
ing parameters. Our purpose here is to relate these vector 
fields with the response coefficients of Eq. (fl5l) . As we will 
show in the next section, using these fields we can con- 
veniently describe the conditions for finite pumping and 
provide a detailed insight into the "stability diagrams" 
for the pumped charge and spin. 

According to Eqs. (fTT)) and (IT8l) . the charge and spin 
pumped in a cycle of the modulation can be written as 
the line integral 



(21) 



Qr = f dx ■ Ar{x) 



For the two-dimensional parameter space studied here, 
spanned by e Xl — (1,0) and e X2 = (0,1), the position 
in the closed trajectory C is indicated by % = X)iXi e i- 
This integral is independent of how fast the path is tra- 
versed, and consequently the total pumped charge (or 
spin) does not depend on the driving frequency as long as 
the adiabaticity condition is fulfilled. The geometric as- 
pects of the problem, entering through the field Ar, are 
certainly of interest However, here they are merely 
convenient auxiliary quantities to analyze the problem 
of interaction-induced pumping. In analogy to classical 
electrodynamics, the vector field 



A R (x) 



E^ 



(22) 



with V = J2i &xi e i> can be interpreted as a pseudo vector 
potential defined in the space of the driving parameters. 
From Eqs. (fTO|) and (|2"0"|) , the components of this vec- 
tor potential are given by the emissivities to the leads, 
i.e. Ai r = dN(r)/dx and Aj r = dS(r)/dx- Therefore, 
the vector potentials describe, respectively, the amount 
of charge and spin entering the leads due to the change in 
the driving parameters. From the form of Ar in Eq. (|22[) , 
we notice that for constant response coefficients, the re- 
sulting pseudo vector potential is just a gauge function 

Ar = J2j ^{ffPj ) which, integrated over a closed tra- 
jectory, gives zero pumping. 

Using Stokes' theorem we write the pumped charge 
and spin in terms of the surface integral 



Qr 



IL 



dS-B R ( X ), 



(23) 



where E is any area in the parameter space encircled by 
C, such that C — <9E. Written in this way, we can imag- 
ine the pumped charge and spin as the flux generated by 
a pseudo magnetic field 

Br(x) = V x Ar( X ) (24a) 
= E V ^?(X) X Vpf(x). (24b) 
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The advantage of this representation is that the pseudo 
magnetic field anticipates the conditions for finite pump- 
ing without referring to the specific details of the model 
and the modulation. The direction of this field is, by 
construction, perpendicular to the plane (i.e. pointing 
outside the parameter space) defined by the driving pa- 
rameters, i.e. Br{x) — &r(x) e xi x e X2- F° r a fixed 
direction of the driving, the sign of the pumped charge 
is given by the sign of Br, which depends on the inter- 
nal details of the system (e.g., Coulomb interaction, cou- 
pling to the leads, spin degeneracy). We stress that this 
interpretation of Qr comes purely from the adiabaticity 
condition in the time-dependent parameters, and should 
not be confused with any other effect due to the external 
magnetic field B eventually present in this setup. 



III. RESULTS 

In this section, we apply the above theory for the spe- 
cific case of a single-level quantum dot slowly driven by 
the gate and bias voltage. In particular, we want to un- 
derstand how the pumped charge and spin are affected by 
the interplay between the local Coulomb interaction and 
nonequilibrium effects induced by a modulation around 
a finite bias. To this end, we introduce expressions for 
the adiabatic charge and spin currents for a general mod- 
ulation and then we take the specific modulation of the 
voltages given in Eq. (@} • 

Our starting point is the description of the time- 
resolved adiabatic charge current given in Eq. (ITU) . This 
is obtained from the explicit calculation of the matrix el- 
ements of the kernels wf and Wj t introduced in the 
previous section. In the evaluation of the response coef- 
ficients related to the occupation of the different states 
of the dot [see Eq. (IT5l) ]. we find out that it is possible to 
combine them in such a way that the adiabatic current 
can be written in terms of the time-derivatives of the in- 



stantaneous average charge (n) 



ft \(i) 



,0) 



pg+pg + 2pga n d 



spin (S'z)j ' = {p^\ — p^' t )/2. We recover the result^ 



where the charge current response coefficients, related to 
the charge and spin in the dot, respectively, read 

(r - 7 ) (r r + 7r ) + p/3 r 



T 2_ J 2 +/3 2 



(26a) 

- r'- 72 ;> — • ( 26b ) 

together with 7 = ^ r 7r, /3 = J2 r Pr, and T = J2 r T r- 
The response coefficients depend parametrically on t 
through the following factors 



7r 



rv 
2 

r, 
2 



E 



T 



f 



U 



?!>[/(£)-/( 



T 

era + U 

T 



(27a) 
(27b) 



Since the adiabatic charge currents flowing into each 
one of the leads are related by particle conservation, i.e. 



y 1 



(a) 



no net charge is accumulated on 



the dot after one period of the modulation. 

Noticeably, the above expression for the adiabatic 
charge current is not restricted to a particular choice of 
the driving parameters. In this sense, Eq. ([25]) is valid for 
arbitrary combinations involving not only the voltages 
V g (t) and Vb(t) but also T r (t), B(t) (in a fixed direc- 
tion), U(t), etc., because the eigenstates of the system 
remain time-independent. For the specific modulation 
we consider here, given by Eq. (J4j) , the time dependence 
exclusively enters in the arguments 



e r<7 (i) = e a (t) - Hr(t), 



(28) 



of the Fermi function f(ui) = [1 +exp(w)] 1 . We now ex- 
tend the above calculation for the adiabatic spin current. 
From the evaluation of the matrix elements of Wj t we 
obtain 



J r,t — fn,t^\ n )t +'?S I i^Wz/t 



(29) 



In this case, the spin current response coefficients are 
1 (r - 7 )& - (T P - j r )/3 



Vs.,t 



2 r 2 - 7 2 + (3 2 

(r + 7 )(r r - 7r ) + /3 r /3 
r 2 - 7 2 + p 2 



(30a) 
(30b) 



For the spin-isotropic quantum dot discussed here, the 
adiabatic spin current only occurs in the presence of an 
external magnetic field. For B — at any time, ro- 
tation symmetry implies (S z )\ 1 ^ — and j3 r — 0, such 
that the current vanishes for all times. The spin cur- 
rents j| a ] and J R a \ are related by spin conservation, i.e. 



y. r Jrt ~ — It (^z)t ' 1 meaning that no accumulation of 
spin is allowed after one cycle of the parameter modula- 
tion. 



A. Displacement current and pseudo gauge 
invariance 

In addition to the tunneling currents we introduced 
above, a displacement current: 



(31) 



generally occurs in the present setup due to the moving 
screening charges <5^. cr in the gate and the leads. These 
arise in response to a variation in the electrostatic poten- 
tial induced by a change of the charge on the dot. To de- 
scribe them, we consider the Coulomb-blockade model^ 
in which the system is represented by the equivalent cir- 
cuit shown in Fig. [5] The screening charges are given by 
the difference 



C (t) = C r [V r (t) - F d ot(Qdot(*))] 



(32) 
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FIG. 2. Schematics of the equivalent circuit for the proposed 
model. Cl , Cr are the capacitances between the dot and the 
leads and C g is the capacitance between the dot and the gate 



electrode, 
lines. 



The quantum dot region is indicated by dashed 



between the applied voltage at the r-lead and the electro- 
static potential inside the dot (in Fig. [21 the region de- 
limited by dashed lines). C r is the capacitance between 
the dot and the r-lead. Due to particle conservation we 
have Qdot(i) = d{h) /dt = — [Ii(t) + Iti{t)\ and the dis- 
placement current reads^ 



vector potential 
A x? . = C r VV r °r'VV r , + , (36) 



which can be imagined as a gauge function, A]<ns = VV>, 
that leaves the total pseudo magnetic field unchanged, 
i.e. Bitot = V x (A Ir + Vip) = V x A Ir . Since the 
pumped charge is due to this pseudo magnetic field, we 
can focus in what follows only on tunneling currents. 

Now that we have the formal expressions for the adi- 
abatic charge and spin currents, we can start with the 
analysis of the pumping for the specific modulation of 
the gate and bias voltage [see Eq. (gj)] and fix the direc- 
tion of the driving by taking <fi g = —ir/2 and 0& = n for 
which the adiabatic dc current is maximal. Motivated by 
a previous study on interaction-induced pumping^ we 
describe in the following sections the pumping of charge 
and spin making use of the vector fields of Eqs. (|22p and 
([Ml) . Since we are interested in the effect of the interac- 
tion U on the pumping, we first consider the case U = 
as a reference. We then analyze the role of the local inter- 
action U > for B = and characterize the mechanism 
of pumping. Finally, in Sec. lIIIDl we extend the descrip- 
tion of the pumped charge and introduce the pumped 
spin by including an external magnetic field B > Q. 



B. Non-interacting quantum dot 



X r dis (i) = G r Vr(t) ~ ^ E [ C r' Vr> («) + U> (<)] , (33) 



where C — Cl + Cr + C g is the total capacitance. As 
in the case of the tunneling currents, we can separate 
the displacement current into instantaneous and adia- 
batic contributions. Since the time-derivatives of the ap- 
plied voltages are of linear order in i7, they do not enter 
in the instantaneous term and 



dis,(i) 



Cr \ ^ T (i) 
' (J r',f 



(34) 



(i) 

However, since j. corresponds to the instantaneous 
steady-state solution [see Eqs. © and ([13"])] . particle con- 
servation implies I L l \ + — and therefore 1^'^ is 
exactly zero for any value of t, as required. For the adi- 
abatic correction to the displacement current we have 



For U = 0, B > we observe that the charge current 
response coefficients [see Eq. (|26p ] are time-independent 
since j r = /3 r = 0. Therefore, the adiabatic response is 
unaffected by the change in the dot occupation during the 
full cycle. This means that the same amount of charge 
that flows into the dot from the r-lead during the loading 
part of the cycle, characterized by d (n)y' /dt > 0, is re- 
turned to the same lead during the unloading. This con- 
stitutes a clear example in which the loading/unloading 
symmetry is preserved after a complete cycle of the driv- 
ing. The analysis is also valid for the adiabatic spin 
current, where the loading/unloading symmetry is mani- 
fested through the average spin inside the dot. Therefore, 
although there is a finite charge and spin current 



t(») _ 

J r,t — 



r,, d 
Ydt 

r dt {bz/t 



(i) 



(37a) 
(37b) 



/d i S ,(a) = Cr y At) _^Y^C r , V r , (t) + 7, 



c 



r(«) 



(35) 



Here one notices that although this current can have fi- 
nite values during the driving cycle, its time-average over 
one period is zero. This becomes evident since particle 



conservation yields / 



(a) T (a) 



L,t 



L R,t 



such that the 



above is a total time-derivative. 

In terms of the vector fields of the previous section, the 
displacement current is related to an irrotational pseudo 



the total pumped charge and spin, obtained after inte- 
grating over the full period, are exactly zero. In terms 
of the vector field introduced in Eq. ((22l) . the constant 
response is then described as an irrotational vector po- 



tential An = V 



with R either 



I r or J ri whose corresponding pseudo magnetic field 
Br = V x Ar : is indeed zero. We emphasize that this 
result depends on the particular choice of the driving pa- 
rameters. Other setups involving time-dependent barri- 
ers yield finite pumping even in the {7 = limit since the 
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response coefficients are time-dependent and the load- 
ing/unloading symmetry is not necessarily preserved*^ 

C. Interacting quantum dot at zero magnetic field 

We now focus on the effect of a finite Coulomb in- 
teraction U 3> T > T by considering B = 0. Once the 
mechanism that generates the pumping is understood, we 
extend, in Sec. IIII D[ the discussion to a finite magnetic 
field. 

In Fig. [3]we show the pseudo magnetic field of Eq. 
in terms of the driving parameters V g and (from hereon 
called 'stability diagram for the pseudo magnetic field for 
the pumped charge'). This map is a convenient tool in 
the characterization of the pumped charge (and later on 
the pumped spin) for different regimes of the applied bias 
and arbitrary asymmetry in the coupling to the leads. 
In this example, we show the pseudo magnetic field Bi L 
related to the adiabatic charge current entering the left 
lead. As observed in Ref. 43, it displays four peaks of lin- 
ear dimension ~ T, located in the vicinity of the meeting 
point of two resonance lines of the differential conduc- 
tance, namely, at the corner points of the regions of stable 
charge in the quantum dot. These lines can be approxi- 
mated by the dot-level resonance conditions fi r = e,e + U , 
respectively, with e = —aV g and r = L,R (grey dashed 
lines in the figure). The low-bias peaks [labeled by (1) 
and (3) in Fig. [5] are dominant, while the high-bias peaks 
[(2) and (4)], emerge only when we include an asymme- 
try A = (r^ — Tn)/r between the two barriers. We first 
investigate the regions in which Bj L is zero and then we 
separate the discussion of the peaks according to the ap- 
plied bias around which the driving takes place. 

Since we consider B — 0, the singly occupied dot states 
are degenerate, i.e. e rcr = e — /i r for a =j, J,, such that 

(i) 

(S z } t =0 and the adiabatic spin current completely 
vanishes. In this case, the adiabatic charge current re- 
duces to 



(a) 



r r + 7 r d 

r + 7 ~dt 



v(0 



W W =2 E p rr/(er/T) 



r + 7 



(38a) 
(38b) 



and admits the following interpretation: As soon as the 
driving passes through a resonance line, the occupation 
in the quantum dot is changed and generates a response 
current flowing from/into the leads. The contribution 
flowing through the r-barrier is then given by the ratio 
(r r + 7r)/(r + 7) between the charge relaxation rate rel- 
ative to the r-lead, i.e. T r + j r , and the total relaxation 
rate r + 7 corresponding to the sum of the two leads. 
When taking the time-average over a complete cycle, the 
occurrence of a dc component of the adiabatic current is 
tied to an asymmetry between the loading and unload- 
ing parts of the cycle. We can find two regimes in which 
this condition is not fulfilled: (i) When the driving is 
far away from any resonance line, the occupation in the 
dot remains constant and the response is exactly zero. 
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FIG. 3. (Color online) Stability diagram for the pumped 
charge: normalized pseudo magnetic field as function of the 
driving parameters V g and V&, for B = 0. Grey dashed 
lines correspond to dot level resonance lines. The chosen 
Coulomb interaction and junction asymmetry are, respec- 
tively, U = 30T and A = 0.25. 



(ii) When the driving crosses a single resonance line, the 
response is symmetric during the loading and unloading. 
Therefore, the same amount of charge that enters the dot 
from the r-lead during the loading returns to the same 
lead in the unloading, such that the dc component of the 
adiabatic current is zero. In terms of the pseudo mag- 
netic field 



(39) 



we observe that far away from any resonance line, the gra- 
dient of the average charge is exponentially suppressed. 
On the other hand, when the trajectory dT, traced by the 
driving parameters only crosses a single resonance line, 
the response coefficient tp% and the average charge (n) ^ 
depend on the same effective parameter and the vectors 
in Eq. ([39]) become parallel to each other. In this sense, 
adiabatic transport along a single resonance line can be 
understood as single parameter pumping, 5 - which is well 
known to give zero contribution to the time-averaged cur- 
rent in the adiabatic limit iSl 

In the remainder of this section, we analyze the dis- 
crete points around which the pseudo magnetic field is 
nonzero and calculate the 2D resonance shape of the re- 
lated pumped charge. To this end, we first notice that 
for the modulation we consider in Eq. (j4]), it is conve- 
nient to write the driving parameters as they enter in 
the arguments of the Fermi function, i.e. 



Xr(t) 



aV g (t) V b (t) 



+ 



(40) 



T 2T ' 

for r = L,R, respectively. Motivated by the location of 
the peaks in the pseudo magnetic field, we investigate 
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first the pumping at low bias, characterized by the peaks 
marked by (1) and (3) in Fig. [3] and then we consider the 
high-bias peaks (2) and (4). Once we know the behavior 
of these peaks, a complete description of the stability 
diagram for the pumped charge can be obtained from 
the symmetries of Bi r with respect to a reflection of the 
applied voltages (see Appendix lAl . 



1. Low-bias regime 

We now calculate the pseudo magnetic field for the 
low-bias peak (1) around the point (aV g ,Vb) = (0,0). In 
this case, we take the gradient of the response coefficient 
and the average charge in Eq. ([39jl with respect to the 
driving parameters of Eq. (|40l) and obtain 



r 



8(r + 7 ) 3 cosh 2 (2^) cosh 2 (2^) 



(41) 



with 7 ~ J2r T rf{Xr) and x = (xl,Xr)- The super- 
script in Bi L labels the origin of coordinates with respect 
to which xl and XR 1& measured. Notice here that the 
sign of the field is independent of the coupling asym- 
metry A. In particular, for the chosen direction of the 
pumping cycle, the positive sign indicates that the re- 
sponse current through the left lead is stronger during 
the unloading part of the cycle. Additionally, the pseudo 
magnetic field decays exponentially for x ^ 1, which, as 
we will show later, implies an asymptotic value for the 
pumped charge when the enclosed area is larger than the 
typical support of the field (~ 5T). 

Now that we know the specific shape of Bi L , we can 
determine the condition at which the pumped charge is 
maximal, i.e. we need to know the position (xl,Xr) 
where the pseudo magnetic field reaches its maximum 
value. The exact position is obtained from the roots of a 
quartic equation (see Appendix[B]), which can be approx- 
imated through the interpolation between the maximum 
points for A = and A = ±1, i.e. 



Xx„R^ln(l+\/3)A 



Ail 



dn 



1 + V33 N 



(1-A 2 ). (42) 



In Fig. U we show the trajectory of the maximum value of 
the field, whose location is given in Eq.[22J as the junction 
asymmetry is swept over the whole range — 1 < A < 1. 
In the symmetric case A = 0, we observe that the peak 
is shifted with respect to the charge degeneracy point 
(aVg,Vb) = (0,0). Its position is given by XL = XR = 
ln[(l + v33)/4] — 0.522 or, in terms of the gate and bias 
voltage: aV g ~ —0.522 T, H = 0. The temperature de- 
pendent shift in the pseudo magnetic field is similar to the 
known shift in the SET peak of the linear conductance 
This last is defined via the instantaneous current by 



G = 



dl 



(0 



dV b 



1-A 2 



1 



v b =o 



8Tl + /(£) cosh 2 (^ 



(43) 
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FIG. 4. Position of the low-bias peak (1) of the pseudo mag- 
netic field as function of the junction asymmetry. The curve 
shown here, given by Eq. (|42l> . cannot be distinguished from 
the exact one. 



In addition to the broadening, the peak in G shows a 
shift Ae = f In 2 ~ 0.347 T which increases linearly with 
T as a consequence of the Coulomb interaction. This 
effect is related to a change in the spin degeneracy of the 
ground state when crossing the charge degeneracy point 
(aV g ,Vb) — (0,0). The pseudo magnetic field can be 
expressed in terms of the linear conductance through 



B\ 



(0,0) 



1 



v b =o 



1 1 - A 2 



T 



32 [1 + /(#)] 3 COsh4 (#)' 



(44a) 
(44b) 



and though the origin of the peak shift in Eq. ([42"T) is the 
same, the above level- dependent prefactor explains the 
different value as compared to the one in G. 

For a finite asymmetry in the tunnel couplings, the 
shift in the gate voltage (see Fig. 0} remains almost 
constant and the peak moves vertically in the stability 
diagram. The extreme cases A = — 1 and A = 1 im- 
ply {aV g ,V b ) ~ (-T/2,T) and (-T/2,-T) respectively, 
such that the peak sits close to the resonance lines, as 
indicated by the arrows in Fig. 2J 

The extension of the above analysis to peak (3) in 
Fig. [3] is straightforward. In this case, we set the ori- 
gin of coordinates at the point (aV g , H) = (U,0) and the 
resulting pseudo magnetic field can then be related to the 
one of Eq. (gTJ| by 



BT' Q \ X L,XR) = -B?f(-XL,-XR)- 



(45) 



The negative sign in the above equation indicates that 
now the loading part dominates the cycle. Due to this 
sign, we now look at the points that minimize the field 
here. According to the inversion of the sign in the driving 
parameters, these points are the ones of Eq. (|4*2"j) but 
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Qff = 1/6 
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FIG. 5. (Color online) Pumped charge in the low bias regime 
as function of the modulation amplitude 8\ f° r several val- 
ues of the junction asymmetry. For A = 0, the regimes of 
small and large driving amplitudes [see Eqs. (|46[) an d (HZ} ] 
are indicated by dashed red and gray lines respectively. 



with opposite sign. Therefore, since we took the origin 
at (aV g ,Vb) = (U,0), the peak (3) is shifted into the 
N = 2 Coulomb diamond, such that for A = 0, the peak 
is located at aV g ~ U + 0.522 T and V b = 0. Finally, 
when taking finite values of A, the shift of the position of 
the peak (3) as a function of the bias is opposite to the 
one observed in the peak (1). In this sense, for A — » — 1 
the peak approaches to (aV g , Vb) = (U + T/2, — T) and 
when A 1 it sits close to (U + T/2, T). 

We now calculate the maximal pumped charge as ob- 
tained from Eqs. (|23"]) - (|2"4]) when the working point is set 
in the position where the peak is maximum [see Eq. (|42[) ]. 
We consider a circular trajectory <9£ for the driving pa- 
rameters as defined in Eq. ([J). The numerical evaluation 

of the pumped charge related to the peak (1) is 

depicted in Fig. [5] as a function of the modulation ampli- 
tude Sx, i.e. the radius of the circle over which the field 
is integrated. We observe that, regardless of the value 
of A or Sx, the sign of the pumped charge is fixed only 
by the direction of the pumping cycle. As we will show 
in the next subsection, this is not the case at high- bias, 
where the pumped charge shows a strong dependence on 
A. 

For small driving amplitudes, i.e. 5% <C 1, the pumped 
charge is proportional to the area E encircled by the 
trajectory dT, of the driving cycle, in agreement with 
Refs. [2l] and HH For an arbitrary junction asymmetry, 
this can be approximated by 

Q?f c(l-X 2 )BTr^ X 2 , (46) 

where B^ 3 * — 0.0106 is the maximum value of the peak 
for A = Q. 

When increasing the amplitude, the driving param- 
eters start exploring regions which are away from the 
crossing point of the two resonance lines and the pumped 
charge no longer follows the above relation. In this large 
amplitude (but still adiabatic) driving regime, charac- 
terized by Sx ^ 1) it is important to remark that the 
adiabaticity condition flSx C T is still preserved, since 



we can always take arbitrary small values for the modula- 
tion frequency without affecting Qj L . When evaluating 
the pseudo magnetic field in Eq. (f3"9"| for Xr ^> 1, the 
exponential decay implies an asymptotic value for the 
pumped charge, as shown in Fig. [5l In this case, the 
particular choice of the working point and the specific 
shape of the trajectory in the parameter space become 
irrelevant as far as the peak is fully contained in S. The 
pumped charge then saturates at^ 



Qi 



(0,0) 



3 1- A 2 
2 9 - A 2 ' 



(47) 



It shows a quadratic dependence on A for A <C 1, in agree- 
ment with Fig. [5] The maximum value of the pumped 
charge, corresponding to A = 0, is 



0(0,0). _l 
W I L I max g j 



(48) 



in units of the electronic charge. To explain this partic- 
ular value and to illustrate the mechanism of pumping 
at low bias we show, in Fig. [51 the time-resolved aver- 
age charge and response coefficient during a cycle of the 
modulation around the working point (aV g , Vb) = (0,0). 
We divide the cycle into four steps corresponding to the 
different regions of the stability diagram visited by the 
driving parameters. Here we consider a symmetric cou- 
pling to the leads (A = 0). The four regions are: 

a) Turning point (meaning that the dot level takes its 
maximum or minimum value) above the chemical po- 
tentials: e > fx r . The 7 r factors are exponentially 
suppressed and ^ — —1/2. 

b) Going down between the chemical potentials: /il < 
e < fiR. Here jl ~ and — T/2, such that ^ = 
-1/3 and ipif = -2/3. 

c) Turning point below: e < fi r . Here 7 r ~ T/2 and 
the same response coefficients are obtained. As in a), 
these arc —1/2. 

d) Going up through /xr < e < the asymmetric 
situation observed in b) is reversed and we obtain 



9,; 



-2/3 and tp}* = -1/3. 



To estimate the pumped charge, we consider the time- 
integral of the current given in Eq. (|38[) noting that, in 
the large amplitude driving regime, the time-derivative 
of (h)i is only nonzero when crossing a resonance line. 
This last derives from the blocking of certain transi- 
tions by the Coulomb interaction. Therefore, we can 
approximate d{h)^ jdt through the difference between 
the asymptotic values obtained at each side of the reso- 
nance, while for the response coefficient we take the av- 
erage between its asymptotic values before and after the 
crossing, i.e. Q { ^ 0) ~ -5/18-5/36 + 7/18 + 7/36= 1/6. 
This corresponds to the maximum possible value of the 
pumped charge when modulating the gate and bias volt- 
age. We observe that for this particular modulation the 
total pumped charge is not quantized. Such quantization 
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us to write 
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FIG. 6. (Color online) Time-resolved average charge (solid 
red) and response coefficient (dashed blue) for A = and 
5% > 1 during a cycle of the modulation. Inset: scheme 
of the explored regions of the stability diagram during the 
pumping cycle around the working point (aV g , Vb) = (0, 0). 



of Qj^ would demand the modulation of an additional 
parameter (e.g. the tunnel barriers) 6 2> and is not what 
we address here. In general, the (measurable) plateau 
value that is reached does, however, provide information 
about the tunnel coupling asymmetry. This is similar to 
the use of noise values in the SET regime^ 

To complete the analysis of pumping at low bias, we 
mention that the pumped charge Q^'°^ related to the 
peak (3) around the point (aV g , Vb) — (C7, 0) reads 



n (u,o) _ _ n (o,o) 



(49) 



where we used the antisymmetric shape of Bi L along the 
particle-hole symmetry point (aV g ,Vb) — (U/2,0) [see 
Eq. (g5j) and Appendix E]. 



2. High-bias regime 

Now we extend the above discussion to the region 
around a large static bias Vb ~ U, i.e. the peak (2) 
in Fig. |31 We evaluate the pseudo magnetic field and 
the resulting pumped charge when encircling the cross- 
ing point [cf. Eq. fl4}] (aV g ,V b ) = (U/2,U), which is 
taken as the new origin of coordinates for the voltages. 
In this regime, the resulting pseudo magnetic field can 
be written in terms of the low-bias field as 



(50) 



This relation between low- and high-bias field is a central 
result of this paper. As a direct consequence of the A 
prefactor, we observe that since |A| < 1, the magnitude 
of the peak (2) is always smaller than the one of peak 
(1) and its sign is uniquely determined by the sign of A. 
For any two modulation curves of the same shape and 
direction, centered around these points and symmetric 
with respect to the \R axis, a change of variables allows 



/ f dSB^ U \x) = \JJ dSB?-°\ X ), (51) 



such that we can calculate QfJ 2 ' U%> 
pumped charge at low bias 



)(U/2,U) 



A i 



? (00)_ 



in terms of the 



(52) 



As noted in Ref. |43|, the mere presence of a pumped 
charge in the high-bias regime indicates an asymmetric 
coupling to the leads. Since for the chosen modulation 



the sign of 



is always positive, the sign of Qj U J 2 '' 



could be used as a quick test to determine which one 
of the two leads is dominating the transport. For a di- 
rect quantitative estimation of A, one simply divides the 
pumped charge at the different bias regimes: 



Tl-Tr Q)' 



((7/2,(7) 



)(0,0) 



(53) 



In particular, this is convenient in the regime of large 
driving amplitudes, where the pumped charge is not af- 
fected by the precise details of the trajectory. 

As compared to the low-bias peak (1), the A prefac- 
tor tells us that the loading/unloading symmetry can no 
longer be broken for symmetric barriers, i.e. A = 0. In 
fact, for this particular case, although there is a change 
in the response coefficient along the pumping cycle, it 
has the same time-dependence as the average charge, i.e. 



1, 



(54) 



and the adiabatic charge current is a total time- 
derivative. Therefore, a finite pseudo magnetic field is 
now a cooperative effect of a change in the response coef- 
ficient and the junction asymmetry: A finite A is required 
in order to have non-parallel gradients in Eq. (|3"9"1) . 

Eq. [50] is also useful in that it allows us to determine 
the maximal pumped charge working point based entirely 
on the low-bias feature. According to Eq. (I5TJ1) . the de- 
pendence on A of this point follows the same condition 
as in the low-bias regime, except for the sign inversion 
of xl and the shifted origin of coordinates. Therefore, 
starting from the position of the peak (1) in Fig. [3J we 
can determine the position of the peak (2) by perform- 
ing first a reflection at the resonance line e = /xj, and 
then a translation by (J7/2, U). Although for symmetric 
junctions (A = 0) there is no peak, its position would be 
shifted in the bias voltage by A 14 = 21n[(l + \/33)/4] T 
with respect to the crossing of the dot-level resonance 
lines at (U/2,U). As soon as we increase A from 0, the 
peak emerges and moves almost horizontally (i.e. along 
the gate voltage axis) towards the e = resonance line 
whereas for negative values of A the peak moves towards 
the e = hl line. 

Finally, the above analysis can be transferred to the 
peak (4) by using the bias voltage symmetry discussed in 
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Appendix [A] In this case, the pseudo magnetic field can 
be written in terms of the low-bias field as 

B?J 2 >- U \xl,X R ) = - ABg^Ott,-**), (55) 
such that the pumped charge at this working point 

Qf L /2 >~ U) = -\Qff\ (56) 

provides an additional and independent quantitative 
measurement of the junction asymmetry. This can be 
used as a cross check on experimental results. 

D. Finite external magnetic field 

We now include a finite external magnetic field B 3> 
T.— In this situation, in addition to the adiabatic charge 



current , a nonzero adiabatic spin current J^" ; also 
flows through the dot. 

In Fig. [7] we show stability diagrams for the pseudo 
magnetic fields for the pumped charge (top panels) and 
the pumped spin (bottom panels) for different junction 
asymmetries. The external magnetic field B = 10 T, 
where T is the thermal energy, now splits the resonances 
of Fig. [3] into further well-separated peaks. The pseudo 
magnetic field is nonzero only around the meeting points 
of two crossing resonance lines. As in the case B = 0, this 
is again because the loading/unloading symmetry is pre- 
served when the driving is far away from any resonance 
line or when it only crosses a single resonance line. A 
simple inspection of Fig. [7J shows that, regardless of the 
value of A, there is no peak in the crossing point at zero 
bias (black arrows in the upper left panel). In Ref. 0, 
the absence of a peak was related to the lifting of the spin 
degeneracy in the states with single occupation. Addi- 
tionally, a strong dependence on A in the high-bias peaks 
[(2) and (3) in Fig. [Jj of Bi L is observed^ In contrast, 
the peaks of B,j L are almost unaffected by A. The dif- 
ference between charge and spin current is particularly 
strong for A = 0, where we find pure spin pumping in the 
high-bias regime: the pumped charge peaks (top cen- 
ter panel) vanish exactly while the ones related to the 
pumped spin (bottom center panel) remain finite. 

To understand the above features of the pseudo mag- 
netic fields and how they affect the pumping, we consider 
first the absence of a peak around (aV g , V&) = (—B/2,0). 
In this regime of the driving parameters, the charge in 
the quantum dot is spin polarized, such that the av- 



O) 



(i) 



and 



erage charge and spin simplify to {n)' t ' = pv- 

(S z )[ ^ = p^\/2, respectively. Therefore, the adiabatic 
currents of Eqs. (f2"5)l and (|2l?|) simplify to 



J 



r(«) 
l L,t 

(a) 
LA 



rid 
r dt 

(a) 



2 J L,t 



(57a) 
(57b) 



Since in this regime there is a single available transition, 
namely, |0) «-> |t), the relative rate at which the dot is 



loaded has to be the same as the one during the unload- 
ing. In the language of the vector fields, the vector po- 
tential A-r associated with these currents is irrotational, 
such that integration over the closed trajectory yields 
zero pumped charge. As compared to the B = result 
at low bias (see Sec. IIII C 1[) . a finite pumped charge re- 
quires not only a modulation encircling the meeting point 
of two resonance lines but also a change in the spin de- 
generacy of the ground state. 

We now investigate the regions in which the pseudo 
magnetic field is nonzero. We consider first the peak la- 
beled by (1) in Fig.JJJ Here, the two transitions |0) «-> \\) 
and |0) «-» 14.) arc enabled by the bias window /i^ — Hr 
such that the spin degeneracy in the N = 1 charge block 
is effectively recovered and the loading/unloading sym- 
metry is again broken. To calculate the pseudo magnetic 
fields associated with the charge (Bi L ) and spin (Bj L ) 
currents, we set the origin of coordinates at the crossing 
point (aV g ,Vb) = (0,5). Although the explicit expres- 
sion for the fields is cumbersome, these show the simple 
relation 



Bt B \x) = 



~<' B) (x), 



(58) 



for arbitrary junction asymmetries. By calculating the 
corresponding pumped charge and spin, and noticing 
that these follow the same relation as the fields, we can 
describe the spin resolved pumped charge through the 
definitions 



Qjt = 
Q T i = 

1 j. 



Qil + 2Q 



Jt 



(59a) 
(59b) 



Since in this regime Qj^ — —Q^ /2, the pumped 
charge is purely given by spin 4- carriers. To estimate the 



(0,5) 



dependence on A, we calculate Q\ L 
amplitude and obtain 

(0,B) = 1 + A 
WlL 3 + A 



for a large driving 



(60) 



Notice that, for the chosen modulation, the sign of the 
pumped charge is always positive, regardless of the par- 
ticular value of A. Since no charge is accumulated after 
one period of the modulation, i.e. + Qi° R B ^ = 0, 

the charge is pumped from the right lead to the left lead. 
In this case, the breaking of the loading/unloading sym- 
metry in both the charge and spin in the dot is mainly 
due to the change in the number of available transitions. 
This affects the value of the response coefficients in such a 
way that the total amount of charge leaving the left lead 
during the loading part of the cycle is smaller than the 
one entering the same lead during the unloading. This is 
illustrated in Fig. HI for the case A = and 5\ ^> 1, where 
we show the average charge and spin in the dot together 
with the corresponding response coefficients. According 
to Eq. (j6"0"|) , the pumped charge for A = corresponds 
to 1/3 of the electronic charge. This limit can be un- 
derstood from the asymptotic values shown in Fig. [51 
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FIG. 7. (Color online) Normalized pseudo magnetic fields for the pumped charge (top panels) and pumped spin (bottom panels) 
for junction asymmetries A = —0.25 (left), A = (center) and A = 0.25 (right). The chosen interaction energy is U = SOT and 
the external magnetic field is B = 10T. 



in a similar way as we did in the previous section for 
The difference now is that the pumped charge is 
not only subjected to a change in the average charge of 
the dot, but also to the variation of the average spin. In 
particular, the loading (unloading) of the spin does not 
necessarily correlate with the loading (unloading) of the 
charge. Therefore we must distinguish the two contri- 
butions. The estimation of the time-derivatives of the 
average charge and spin, together with the average value 
of the response coefficients each time the driving param- 
eters cross a resonance line (gray dashed lines in Fig. [5]) 
yields 



q(o,b) = 1 / 48 + 15 /48 = i/ 3) 



(61) 



with the first term related to the variation of the average 
charge and the second to the average spin. 

The above discussed value 1/3 for the maximal 
pumped charge can be increased or decreased depending 
on the sign of the junction asymmetry [see Eq. (pO)) ]. For 
positive A, the maximal pumped charge increases from 
1/3 to 1/2. However, in the extreme case A ~ 1 the (neg- 
ative) peak (2) approaches the region of integration and 
its contribution can no longer be disregarded. 

We now consider the pumped charge and spin in the 
high-bias limit, characterized by the peaks (2) and (3) 
in Fig. [7J To keep the notation simple, we now use the 
index (2) to indicate that the origin of coordinates is set 
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FIG. 8. (Color online) Time-resolved average charge (solid 
black) and spin (solid blue) together with the charge current 
response coefficients cp 1 ^ (dashed red) and tp^ t (dashed green) 
for A = and 8\ S> 1 during a cycle of the modulation around 
the working point (0, B). 



at the point (aV g , Vb) — ((U — B)/2, U). In the regime 
of large driving amplitude we obtain 



(2) 1 (1 + A) 2 (3-A) 
■ Jl 2 (3 + A)(3 + A 2 ) : 



-2XQf. 



(62a) 
(62b) 
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Therefore, since the pumped spin is always positive, the 
resulting pumped charge changes its sign in the symme- 
try point A = 0. Notice that the regime of validity of 
the above limit includes |A| < 1. For A ~ 1 we should 
consider the contribution from the peak (1), such that 
both the pumped charge and spin go to zero. By using 
Eq. (|59[) . the spin resolved current shows to be decom- 
posed into contributions from spin f carriers flowing from 
the right lead to the left lead and spin J, carriers flowing 
in the opposite direction. The ratio between these two 
contributions depends on A via 



1 i 



■> j -V 



1- A 
1 + A' 



(63) 



such that for A < (> 0) the transport is dominated by 
spin f (J,) carriers. Remarkably, for a symmetrically cou- 
pled dot we obtain pure spin pumping, since the two con- 
tributions are exactly opposite and therefore the pumped 
charge is zero. In this case, the limit S\ 3> 1 yields 



Qf, = 



(64) 



which can again be understood in terms of the asymptotic 
values of the response coefficients and the instantaneous 
average charge and spin as we previously demonstrated 
for Q ( £ B) [see Eq. flST}]. 

Finally, we analyze pumping for large driving ampli- 
tudes around the peak (3), given by the crossing of the 
resonance lines at (aV g) H) = (U/2, U + B). Now the re- 
lation between the pumped charge and spin is non-trivial: 



1-A 2 
3 + A 2 ' 
1 3 + 4^ 

4 V 



gg> = a 



^ = 7(1 -A 2 ) 



(3 + A 2 ) 2 



(65a) 
(65b) 



It implies a positive sign for the pumped spin while the 
pumped charge again changes its sign with the sign of A. 
Despite this, the spin resolved pumped charge, calculated 
by Eq. obeys Q ( V{X) = -Q ( 5(-A). As in the pre- 

viously studied regime around the peak (2), symmetric 
coupling to the leads yields a pure pumped spin, which 
in this case is 



12 



(66) 



IV. CONCLUSIONS 

We investigated adiabatic charge and spin pumping 
through an interacting quantum dot driven out of equilib- 
rium by a non-linear bias voltage and time-dependent pa- 
rameter modulations. We showed that, regardless of the 
specific modulation, the time-resolved adiabatic charge 
and spin currents can be interpreted as the response to 
a perturbation in the instantaneous average charge and 
spin due to the variation of the driving parameters. This 
allowed us to identify the charge and spin emissivities 
to the leads in the context of interacting systems in the 
non-linear transport regime. 

For the specific case of a modulation of the gate and 
bias voltage, we discussed the conditions for interaction- 
induced pumping in terms of the properties of vector 
fields associated to the adiabatic charge and spin cur- 
rents. We observed that the Coulomb interaction is 
crucial since it gives the rotational contribution to the 
pseudo vector potential that cannot be gauged away: It 
generates a nonzero pseudo magnetic field and in conse- 
quence a finite pumped charge. For a single-level quan- 
tum dot, we explored the stability diagram associated to 
this vector field for arbitrary bias and junction asymme- 
try. The shape of the pseudo magnetic field reflects the 
two-parameter condition required for adiabatic pumping, 
such that it shows a maximum whenever two lines of the 
usual dl/dV stability diagram meet. The analytic ex- 
pressions for the pseudo magnetic field and the pumped 
charge enable detailed fitting of experimental results. For 
low bias voltages, the pumping mechanism is dominated 
by the change in the response coefficient when exploring 
different regions of the pumping cycle. In constrast, in 
the high-bias regime, the finite pumped charge is gener- 
ated by the cooperative effect of the above mechanism 
and the junction asymmetry. This allows for a direct 
quantitative determination of the junction asymmetry by 
two measurements of the pumped charge. 

The role of the external magnetic field was found 
to be twofold: First, it restores at low bias the load- 
ing/unloading symmetry previously broken by the lo- 
cal interaction in combination with the spin-degeneracy. 
This is evidenced by a suppression of the pumped charge 
at zero bias. Second, in addition to the pumped charge, 
a pumped spin arises once the spin degeneracy is effec- 
tively recovered through an applied bias. In particular, 
the weak dependence of the pumped spin on the junction 
asymmetry allows for pure spin pumping in the high-bias 
regime. 



The difference with the situation around the peak (2) is 
that now, for A < (> 0), the transport is dominated by 
I (t) carriers. 

Finally, we note that the above discussion of peaks 
(l)-(3) can be extended to all remaining peaks in Fig. [7] 
by using general gate- and bias- voltage symmetries of the 
problem in an external magnetic field which are presented 
in Appendix [Al 



ACKNOWLEDGMENTS 

We acknowledge helpful discussions with M. Biittiker, 
F. Haupt, R.-P. Riwar and L. E. F. Foa Torres, and finan- 
cial support from the DFG under Contract No. SPP-1243 
and the Ministry of Innovation NRW. 



15 



Appendix A: Applied voltage symmetries 

In order to completely characterize the adiabatic trans- 
port of charge and spin along the full stability diagram, 
we derive the reflection symmetries of the pseudo mag- 
netic fields Bi L and Bj L . Specifically, we study the be- 
havior of Br with respect to: (i) the reversal of the bias 
voltage Vb — > — Vb and, (ii) the reversal of the gate volt- 
age aV g — > U/2 — aV g . To simplify the notation, we use 
v = Vb/2. These symmetries derive from the interchange 
of the two electrodes and the particle-hole symmetry of a 
single interacting level. Although not exact anymore for 
multi-level quantum-dots, we expect similar qualitative 
correspondences to hold in general. 



1. Bias voltage reversal 

We write the coupling strength in terms of the junction 
asymmetry by r r (A) = T(l + a r \)/2, with a r — ± for 
r = L,R respectively. The factors in Eq. ([27)1 then obey 
the following relations 



7 r (e,-u,A) = Jr(e,v, -A), 
j3 r (e,-v,X) = 0r(t,v,-\), 



(Ala) 
(Alb) 



where f = R(L) for r = L(R) and e = — aV g . By plug- 
ging these into the definitions of the charge current re- 
sponse coefficients [see Eq. (|26|) ] we obtain 



<p%- (e, -v, A) = -l-tp%-(e,v,-\), 
</?£:(£, --M) = -^(e,u,-A), 



(A2a) 
(A2b) 



and for the spin current response coefficients of Eq. ([50)1 
we have 

<pfr(e,-v,\) = -<p£(e,v,-\) i (A3a) 
p£ (e, -v, A) = — 1 — p£ (e, v, -A). (A3b) 

Now we repeat this analysis for the average charge and 
spin. The explicit expressions of these two quantities for 
arbitrary U and B can be calculated from the instanta- 

(i) 

neous occupation probabilities pj obtained as the solu- 
tion to Eq. ©. These write as follows 



(S, 



E rg r r /(£ rff )(r-2 7(T ) 
r 2 - 7 2 + [i 2 

2 



r 2 



7 2 + r 



(A4a) 
(A4b) 



where a — —a, f (lu) = /(— oj), with f(uj) = [1 + 
exp^/T)]- 1 and la = £ r T r /2[f(e ra ) - f(e ra + U)}. 
For these averages, a change in the sign of the bias volt- 
age is equivalent to an inversion of the tunnel barriers, 
i.e. 



(n)» (e,-v,X) = (n) W (e,v,-X) 
(S z f ] (e,-v,X) = (S z f ) (e,v,-X). 



(i) 



(A5a) 
(A5b) 



In the calculation of the gradients of such quantities, we 
notice that 

(d e <p) (e, -v, A) = - (d € <p) (e, v, -A), (A6a) 

(d v cp) (e, -v, A) = (d v ip) (e, v, -A), (A6b) 

for the response coefficients and 

(d e (R) il) ) (e, -v, A) = (d e (R) (l) ) (e,v,-X), (A7a) 

(d v (R) il) ) (e, -v, A) = - (d v (R) il) ) (e, v, -A), (A7b) 



for the average charge and spin. The pseudo magnetic 
fields then write as follows 

B lL (e,-v,X)=B lL (e,v,-X), (A8a) 
Bj L (e,-v,X) = Bj L (e,v,~X). (A8b) 

These two relations can be directly checked by compar- 
ing the left and right panels of Fig. [7J For A = (see 
central panels) the above equations imply a symmetric 
shape of the pseudo magnetic fields around the zero-bias 
axis. 



2. Gate voltage reversal 

We calculate now the pseudo magnetic fields for an 
inversion of the gate voltage around V g = U/2. In this 
case, we take as the new origin of coordinates the point 
(aV g ,V b ) = (172,0) such that the factors in Eq. ([27|) 
write 



7,.(e, v, A) 



E£# ) 



[/(« 



_ u 

ra 2 



)-/(<** + %)] 



(A9a) 



&(e, v,X)=Yl - f ) - /(^ + £)] , 

(A9b) 

where e ra = e — a r v — aB/2. For these new coordinates, 
the inversion of the gate voltage is then given by e — > —e, 
and the above factors are transformed according to 



j r (-e,v,X) = 7 f (e, v, -A), 
f3 r (-e,v,X) = -p f (e,v,-\). 



(AlOa) 
(AlOb) 



In consequence, the response coefficients obey the follow- 
ing relations 

(-e, v, A) = -1 - y% (e, v, -A), (Alia) 

y& (-e,v, A) = p£(e,v,-A), (Allb) 

^(-e,u,A) = ^(e,«,-A), (Allc) 

p£ (-e, u, A) = -1 - p£ (e, v, -A). (Alld) 
For the average charge and spin we obtain 

(n) (l) (-e,v,X) = 2-(nf ) (e,v,-X), (Al2a) 

(§ z f\-e,v,\)={§,)®(e,v,-\), (A12b) 
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such that the pseudo magnetic fields present different 
symmetries with respect to a change in the gate voltage, 



Bj L (-e,v,X) = —Bi L (e,v,—X), 
Appendix B: Maximum values of Bi L 



(A13a) 
(A13b) 



By solving the first equation, we obtain 



V = 



(1 - X)x(x - 2) 
2a; 2 -2(3 + A)a; + l + A' 



(B3) 



and plugging this into the second equation we obtain the 
following quartic equation for x: 



In this section we calculate the point in which the 
pseudo magnetic field is maximum, i.e. the position of 
the peak (1) in Fig. [3] Our starting point is the explicit 
form of the pseudo magnetic field at low bias. According 
to Eq. (|4ip. the condition for an extremum point (xl, Xr) 
of the pseudo magnetic field is the following: 



tanh 



(f 



-3(l + X)d x J L 



tanh(f) 



2+(l + \)f L + (l-\)fn 

~3(l-\)d XR f R 
2 + (1 + A)/ L + (1 - X)f L ' 



(Bla) 
(Bib) 



where f r = f(xr), r = L,R. Now we take the replace- 
ments x = 1 + e XL and y = 1 + e XR , such that x, y > 1 
and the above equations read 



3(1 + % 
3(1 - X)x 



x-1 

y-2 



y 



l 



(2xy + x + y + X(y - x)), (B2a) 
(2xy + x + y + X(y - x)). (B2b) 



ax + bx 3 + ex 2 + dx + e = 0, 



with 



a = 2A - 6, 
z = 2A 2 - 34, 
(A + l) 2 . 



c 



b = 
d = 



-A 2 -4A + 29, 
-2(A+1) 2 , 



(B4) 



(B5a) 
(B5b) 
(B5c) 



For the particular case of symmetric junctions, the so- 
lutions compatible to the condition x > 1 are x = 
(5 + \/33)/4 and x = (7 + \/37)/6 respectively. Now 
if we use these values in Eq. (IB3|) , the only solution that 
fulfills the condition y > 1 is y = x = (5 + v33)/4. 
Therefore, for symmetric junctions, the maximum of the 
pseudo magnetic field is located at the point xl = XR = 
ln[(l + V33)/4]. 
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